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ABSTRACT 



A summary of Bragg deflection theory and various approaches 
to direct detection acoustooptic spectrum analysis (AOSA) 
modeling is presented. A suitable model is chosen and 
extended to include the effects of diffraction efficiency, 
transducer efficiency, irradiance profiles of incident laser 
illumination, aperture size of the Bragg cell, and the acoustic 
attenuation experienced by the acoustic wavetrain generated by 
the input r-f signal. A Fortran program is developed to model 
the AOSA and predict the output image plane intensity profiles. 
A second version of the program includes a time variable per- 
mitting dynamic simulation of the system response. 



4 



TABLE OF CONTENTS 



I. INTRODUCTION 11 

A. OVERVIEW 11 

B. THESIS OBJECTIVE 12 

C. ACOUSTOOPTICAL SPECTRAL ANALYSIS 13 

D. SPECIFIC GOALS 14 

1. Bragg Diffraction 14 

2. Computer Model 17 

E. PERSPECTIVE 18 

II. ACOUSTOOPTICAL SPECTRAL ANALYSIS (AOSA) MODELING 20 

A. ACOUSTOOPTIC INTERACTION THEORY 20 

B. ACOUSTOOPTIC SPECTRUM ANALYSIS THEORY 26 

C. AOSA MODEL CONSTRUCTION 32 

1. Irradiance Profiles 32 

2. Aperture Size 39 

3. Acoustic Attenuation 44 

4. Diffraction Efficiency 4 6 

5. Transducer Efficiency 57 

D. AOSA MODEL UTILIZATION AND EXTENSION 59 

III. SUMMARY AND CONCLUSIONS 80 

A. SUMMARY 80 

B. AREAS FOR FURTHER STUDY 80 

1. Nonlinearities • - 8 0 

2. Time Varying Signals 81 

3. System Components 81 



5 



C. CONCLUSIONS 



82 



APPENDIX A Acoustooptic Spectrum Analyzer Program 8 3 

APPENDIX B Dynamic Acoustooptic Spectrum Analyzer 

Program 99 

APPENDIX C Arbitrary Profile Solution Details 108 

LIST OF REFERENCES 110 

INITIAL DISTRIBUTION LIST 113 



6 



LIST OF TABLES 



I. Beam Profiles F(z) and Their Corresponding Angular 

Spectral Functions G(B q ) 25 

II. AOSA Modeling Parameters 60 



7 



LIST OF FIGURES 



1. Direct Detection Acoustooptical Spectrum Analyzer 15 

2. Oscilloscope Display of a 70 MHz Monotonic Input 

Signal to the NPS Acoustooptic Spectrum Analyzer 16 

3. Bragg Diffraction Geometry 22 

4. Bragg Cell Operation 27 

5. Classical Ray Trace Diagram for an Acoustooptical 

Spectrum Analyzer 29 

6. Optical Plane Schematic for an Acoustooptical 

Spectrum Analyzer 29 

7. Block Process Element Diagram for an Acoustooptical 

Spectrum Analyzer 30 

8. Bragg Deflection Beams for Various Incident Irra- 

diance Profiles 34 

9. Comparative Image Plane Intensity Profiles (Magni- 
tude) for Gaussian and Plane Incident Beams 37 

10. Image Plane Intensity Profile (Magnitude) Produced 

by a Monotonic 70 MHz Signal with an Incident 
Gaussian Beam and a Simple Transducer 38 

11. Image Plane Intensity Profile (Magnitude) Produced 

by a Monotonic 70 MHz Signal with an Incident Plane 
Wave and a Simple Transducer 4 0 

12. Attenuation Factor as a Function of Frequency for 

Lithium Niobate Transducer Deflector 45 

13. Normalized Transmission as a Function of Frequency 

for Lithium Niobate Transducer Deflector 47 

14. Diffraction Efficiency Bandshape 4 9 

15. First-order Stepped Array Transducer 52 

16. AOD-70 Bragg Deflector Efficiency Curve 54 

17. AOSA Implemented Diffraction Efficiency Bandshape 55 



8 



18. Image Plane Intensity Profile (Magnitude) Produced 

by a 5 MHz Incremental Comb Signal with an Incident 
Plane Wave and a Simple Transducer 56 

19. Image Plane Intensity Profile (Magnitude) Produced 

by a 5 MHz Incremental Comb Signal with an Incident 
Plane Wave and a Stepped Array Transducer 58 

20. Input RF Signal Magnitude for a 70 MHz Monotone 63 

21. Image Plane Intensity Profile (Decibel) Produced by a 
Monotonic 70 MHz Signal with an Incident Gaussian 

Beam and a Stepped Array Transducer 64 

22. Image Plane Intensity Profiel (Decibel) Produced by a 
Mono tonic 70 MHz Signal with an Incident Plane Wave 

and a Stepped Array Transducer 65 

23. Image Plane Intensity Profile (Magnitude) Produced by 

a Monotonic 70 MHz Signal with an Incident Plane Wave 
and a Stepped Array Transducer (Dynamic Model) 56 

24. Image Plane Intensity Profile (Magnitude) Produced by 
a Monotonic 70 MHz Signal with an Incident Gaussian 

Beam and a Stepped Array Transducer (Dynamic Model) 58 

25. Input RF Magntidue for a FM Modulation Signal with 

6 = 8.0, f = 5 MHz 59 

m 

26. Image Plane Intensity Profile (Magnitude) Produced by 

a FM Modulation Signal, 8 = 8.0, f =5 MHz, with an 
Incident Plane Wave and a Stepped Xrray Transducer 70 

27. Image Plane Intensity Profile (Magnitude) Produced by 
a FM Modulation Signal, 8 = 8.0, f =5 MHz, with an 
Incident Gaussian Beam and a Stepped Array Transducer- 71 

28. Image Plane Intensity Profile (Magnitude) Produced 

by a FM Modulation Signal, 8 = 8.0, f = 5 MHz, with 
an Incident Plane Wave and a Stepped Array Trans- 
ducer (Dynamic Model) 72 

29. Image Plane Intensity Profile (Magnitude) Produced 

by a FM Modulation Signal, 8 = 8.0, f = 5 MHz, with 
an Incident Gaussian Beam and a Stepped Array 
Transducer (Dynamic Model) 73 

30. Image Plane Intensity Profile (Magnitude) Produced 
by a Linear Chirp with Af =80 MHz, f =50 MHz, 

t = 2 fisec with an Incident Plane Wave 0 and a Stepped 
Array Transducer 74 



9 



31. 



Image Plane Intensity Profile (Magnitude) Produced 
by a Linear Chirp with Af =80 MHz, f =50 MHz, 
t = 2 ysec with an Incident Gaussian Beam and a 
Stepped Array Transducer 75 

32. Image Plane Intensity Profile (Magnitude) Produced by 

a Linear Chirp with Af m = 80 MHz, f =50 MHz, t = 2 
ysec, with an Incident Plane Wave and a Stepped Array 
Transducer (Dynamic Model) 77 

33. Image Plane Intensity Profile (Magnitude) Produced by 
a Linear Chirp with Af = 80 MHz, f =50 MHz, 

t = 2 ysec with an Incident Gaussian Beam and a 

Stepped Array Transducer (Dynamic Model) 78 

34. Image Plane Intensity Profile (Decibel) Produced by 

a Monotonic 70 MHz Signal with an Incident Plane Wave 
and a Stepped Array Transducer (Dynamic Model) 79 



10 



I. 



INTRODUCTION 



A. OVERVIEW 

Assessment of emergent technologies leads to a continuing 
focus on acoustooptical techniques for spectral analysis and 
filtering/excision [Ref. 1] . Integrated optics with distinc- 
tive advantages in size and weight easily meet the require- 
ments for radar warning receivers. Bulk wave devices offer 
superior dynamic range, bandwidth, and resolution over inte- 
grated optical devices. Foreseeable applications include 
surveillance and electronic intelligence. 

An acoustooptical spectrum analyzer transforms a wideband 
time domain signal into an optical spatially detected fre- 
quency domain signal. Removal of narrowband interference 
may be accomplished by optical excision techniques, to selec- 
tively notch filter unwanted spectral components. U. S. 
military applications in this area could greatly benefit 
from the speed of response and high quality interference re- 
jection capabilities possible with these techniques. The 
potential of these techniques was recognized and led to the 
establishment of a joint NAVELEX 350/DARPA program to examine 
these areas. 

The technical objectives of this effort were: 

1. Survey potential applications and analyze technical 
requirements imposed by each. 

2. Analyze alternative optical excision configurations 
and components. 
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3. Develop and procure major system components. 

4. Fabricate and measure alternative broadband excisor 
configurations to verify analytical results. 

5. Fabricate and field test a brassboard model. 

6. Design a preliminary advanced development model. 

The principal civilian contractor for this project is PROBE 
Systems, Inc. PROBE, assisted by the developmental efforts 
at the Naval Research Laboratory (NRL) and Naval Postgraduate 
School (NPS ) , was tasked with the attainment of these objec- 
tives. NRL's specific program assistance was the development 
of a photodichroic optical clipper to be evaluated by PROBE. 

The research summarized in this thesis is a continuation of 
research initiated as part of the NPS contribution to the 
development program. 

B. THESIS OBJECTIVE 

The objective of this work effort was the attainment of 
an understanding of the principles of direct detection acousto- 
optical spectrum analysis. Previous work at the Naval 
Postgraduate School resulted in the construction by Regan 
[Ref. 2] of a working laboratory acoustooptical spectrum 
analyzer and the subsequent modification of the existing 
equipment by Smith [Ref. 3] to function as an acoustooptical 
spectrum excisor (AOSE) . The acoustooptical spectral excisor 
utilized optical excision, the effect of an object placed 
in that part of the optical path where the undesirable 
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frequency components have been spatially spread. The optical 
excision research included propagation considerations and 
diffraction by the obstacle. The effectiveness of the optical 
excisor in removing interference is dependent upon these 
factors and more fundamentally upon the ability of the de- 
flected beam irradiance pattern to accurately represent the 
input r-f signal spectral components. 

C. ACOUSTOOPTICAL SPECTRAL ANALYSIS 

An acoustooptical spectrum analyzer (AOSA) consists of 
an acoustooptical light deflector (a Bragg cell) , a coherent 
light source (a laser) , and a detector. An input r-f signal 
excites a piezo-electric transducer mounted on the acoustic 
medium on the Bragg cell. This launches an acoustic wave- 
train the length of the cell which modulates the refractive 
index of the interaction medium by the elasto-acoustic 
process. The modulated interaction material functions as 
a moving diffraction grating for the incident laser light. 

A one-dimensional angular pattern of diffracted beams is 
produced and focused on a photodetector array. The pattern 
replicates the spectral content of the original r-f signal. 

The intensity of the individual beams is directly proportional 
to the r-f signal component intensity. The deflection 
angles of the individual beams are proportional to the r-f 
signal frequency component. The detector array photosites 
video detect the incident beam patterns and, in a square law 



13 



fashion, convert the spectral information for electronic 
transmission to a display oscilloscope. Figure 1 shows the 
NPS laboratory AOSA. Figure 2 shows a typical oscilloscope 
display output from the system. Korpel [Ref. 4] includes a 
basics discussion on the theory of Bragg cell applications 
for acoustooptical spectrum analysis in an overview on 
acoustooptics . 

D. SPECIFIC GOALS 

There were three initial goals for this work. The first 
was the study of the factors affecting the Bragg cell 
deflection of incident laser light. The second was the 
analysis of the effect of these factors on the replication 
of the input r-f signal spectrum in the deflected beam 
pattern. The third was the development of a Fortran program 
that would model the acoustooptic spectrum analysis process. 

1. Bragg Diffraction 

The study and analysis of the factors affecting 
Bragg deflection was focused towards an understanding of the 
ability of the deflected beam pattern to faithfully replicate 
the spectral components of the input r-f signal. This 
knowledge was considered essential in order to construct a 
viable computer simulation model. A study of Bragg diffrac- 
tion and the various approaches to modeling of the deflected 
beam irradiance profiles was undertaken. Previous works at 
NPS by Regan [Ref. 2] and Smith [Ref. 3] were used as 
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Figure 2 




Oscilloscope Display of a 70 MHz Monotonic Input 
Signal to the NPS Acoustooptic Spectrum Analyzer 



16 




starting points. Bragg diffraction is traditionally theo- 
rized by use of complex coupled mode wave equations. The 
complexity of these equations rendered them unsuitable for 
the first order simulation sought. A study of other 
approaches was conducted to arrive at a method that would 
provide both adequate accuracy and a reasonable level of 
complexity. 

2 . Computer Model 

An approach to the modeling of the AOSA was selected 
based on the work of Hecht [Ref. 5] . He shows the image 
plane intensity as the squared magnitude of the convolution 
of the spatial Fourier transform of an illumination function 
and the transform of the converted acoustic wave. Some 
modifications in the form of transmission factors were applied 
by algebraic multiplication. Hecht' s illumination factor is 
a spatial weighting function composed of the optical ampli- 
tude profile, the optical aperture size of the Bragg cell, 
and the acoustic attenuation of the interaction material. 

The added transmission factors were the diffraction efficiency 
and the transducer efficiency. 

A Fortran program was developed to model the AOSA 
given the specific physical parameters of the system and 
desired simulation requirements such as laser wavelength and 
laser beam profile. This program was used to study the 
quality of the replication of the r-f signal spectral components. 
The major thrust of this effort was to produce this model for 
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utilization with previously developed optical excision 
programs. While the program is configured to simulate 
various component deflectors, it was structured about the 
AOD-70 model Bragg deflector, manufactured by the 
IntraAction Corp. , which is scheduled for future installa- 
tion in the NPS laboratory AOSA. The AOD-70 Bragg deflector 
is constructed with high quality flint glass and Indium 
bonded Lithium Niobate stepped array transducers. It has 
a 50 MHz to 90 MHz bandwidth and a center frequency of 70 MHz. 

During the development of the program it became 
apparent that the addition of a third dimension, time, to 
the model would significantly enhance the ability of the 
model to demonstrate the response of the AOSA under various 
simulation conditions. A second Fortran program was written, 
essentially a variation of the first, to provide this dynamic 
capability. 

E. PERSPECTIVE 

The application of both AOSA and AOSE techniques in wide- 
band, high density signal environments is dependent on the 
accuracy of the diffracted beam profile in its replication 
of the input signal spectra. The classical model for 
determination of the diffracted beam profile is presented 
by Hecht [Ref. 5] . In a simplified form the diffracted beam 
profile is the convolution of the spatial Fourier transform 
of the illumination field incident on the Bragg cell and the 
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spatial transform of the generated acoustic wave within the 
cell. In detail both the illumination function and the 
acoustic wave are subject to modifying factors. The signifi- 
cant factors include diffraction efficiency, transducer 
efficiency, the irradiance profile of the incident laser 
light, the aperture size of the Bragg cell, and the acoustic 
attenuation experienced by the acoustic wave as it transits 
the cell. This thesis is an attempt to provide a more 
detailed model by incorporating these factors in a first 
order simulation to aid in the study of acoustooptic spectrum 
and excision techniques. 
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II. ACOUSTOOPTICAL SPECTRAL ANALYSIS MODELING 



A. ACOUSTOOPTICAL INTERACTION THEORY 

As an acoustic wave propagates in an isotropic material, 
it produces periodic modulation of the index of refraction 
by means of the elasto-acoustic effect. The interaction of 
incident light waves with acoustic waves is termed acousto- 
optical coupling and is traditionally modeled utilizing 
coupled wave equations. The basic theory of acoustooptical 
interaction in isotropic media was well developed prior to 
the advent of the laser. An excellent discussion of early 
theoretical work can be found in Ref. 6 by Born and Wolf. 

The advent of the laser stimulated research in acoustooptic 
interactions as devices were developed for the deflection 
and modulation of incident laser light. Development of 
superior acoustooptic materials and highly efficient broad- 
band transducers have given impetus to rapid advances in 
these acoustooptic applications. 

Various models describing acoustooptic interaction have 
been proposed. A thin grating model (as it assumes a negli- 
gible depth for the sound field) was proposed by Rosenthal 
[Ref. 7], Korpel and associates [Ref. 8], and Lambert [Ref. 9], 
The thin grating model is essentially a first order solution 
of the diffracted light problem. A second approach known as 
the physical optics model incorporates the Green's function 



20 



integral, as proposed by Gordon [Ref. 10]. An extremely 
rigorous model, the mathematics involved are quite complex. 

A third approach is termed the phenomenological geometrical 
optics model as proposed by Johnson [Ref. 11] and ranks midway 
in complexity between the thin grating and physical optics 
models . 

For the purposes of this work a highly satisfactory 
approach to acoustooptical diffraction is found in a series 
of works by Chu, Kong, and Tamir [Refs. 12-16]. They have 
considered diffraction of arbitrary profiled optical beams 
by periodically modulated layers and have achieved closed 
form solutions for the resulting diffracted far-field 
intensity patterns. A summary of their results is presented 
below. 

As shown in Fig. 3, a bounded beam is incident on a 
periodically modulated dielectric layer with permittivity of 
the form 

e(z) = £ 2(1 + M cos 2irz/d) (1) 

where 

£2 = the relative permittivity of the slab in the 
absence of modulation 
M = the modulation index 
d = periodicity 

The slab has thickness L and is bounded by a dielectric 
medium with relative permittivity £1 for x^O and £3 for x^L. 
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Figure 3. Bragg Diffraction Geometry (after Ref. 12) 
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The electric field of an incident beam may be represented as 



E inc (x ' z) = f G(6 o } exp [ j(C i 1] 



x + B q z) 



] 



( 2 ) 



with 



B q = (2it/X) \/i7 sin 0 



(3) 



?i 1) = (2Tt/X)Ve7 cos 0 



= ^ (2tt/X) 2 £i - B 0 : 



(4) 



G(B q ) is the angular spectral amplitude of the incident beam 
at the entrance plane x=0 . 



G(8 o } = 2? / E inc (0 ' z) exp (_j 6 o z) dz 



= 7{E inc (0 < z)} 



(5) 



Letting 



E inc ^°' z) = F < z ) exp (jbz) 



( 6 ) 



with 



b = (2 tt/X) ViT sin 9 



( 1 ) 



= (2ir/X) VTI sin 9 



(3) 



(7) 



where 



9^^ = angle of incidence of beam axis 
(3) 

9^ = angle of refracted beam 



F (z ) = beam-profile at x=0 
then Equation (5) becomes 
G ( So) = 



/ 



F (z ) exp Oj (S Q -b) z]dz 



( 8 ) 
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Various beam profiles F(z) and their corresponding angular 

spectral functions G(6 q ) are listed in Table I. 

For an optical beam at exactly Bragg incidence we have 

e (1) = 8 b (1) and 9< 3) = e< 3 », and 
o B o B 

b = (2 tt/X)\/£T sin 9 ^ 

B 

= (27T/X)VeT sin 

= ir/d (9) 

The far field pattern for the Bragg deflected beam is given 
by 

P ( 0 ) = G(4>) T(<|>) (10) 

where 

4> = ( 2 t r/X) VeT sin 9 + TT/d 

T(<j>) = is the complex transmission coefficient. 

A description of this coefficient may 
be found in Appendix C. 

Thus the far-field pattern for a Bragg deflected beam is 
essentially the Fourier transform of the aperture field at 
the boundary of the layer modified slightly by a transmission 
factor. 

While the results obtained by Chu, Kong, and Tamir are 
extremely useful in approaching the construction of a model 
for the direct detection acoustooptical spectrum analyzer, 
they are predicated on the periodic modulation of the layer 
by a single frequency acoustic wave. 
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Beam Profiles F(z) and Their Corresponding 
Angular Spectral Functions G(6 ) [after Ref. 
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B. ACOUSTOOPTIC SPECTRUM ANALYSIS THEORY 



The acoustooptic spectrum analyzer consists of only four 
components; a light source, Bragg cell, transform lens, and 
a photodetector array. The basis for acoustooptic spectrum 
analysis is the frequency dispersion of acoustooptic diffrac- 
tion. Essentially each frequency component of the input 
signal will cause a deflection of the incident light beam. 

A typical interaction geometry is shown in Fig. 4 [Ref. 17]. 

A summary of the description of this interaction geometry 
from Ref. 17 is presented. 

Two discrete frequencies are shown with corresponding 
acoustic wavelengths A i and A 2 . These propagate in the x 
direction with velocity V . The medium is optically iso- 
tropic with an index of refraction n. The thickness of the 
ultrasonic grating along the x axis is L. A collimated 
beam of light of freespace wavelength A is incident on the 
sound field at an angle 9^ from the z axis in the xy plane. 

A sound wave at frequency f interacting with the incident 
light wave will generate a principal diffracted beam separated 
from the incident beam by an angle 0^ + 0^ by the well known 
Bragg law; 
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Figure 4. Bragg Cell Operation 
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where 9 . + 6 , < 0.1 rad 

l a 

External to the medium 




~ 2 sin 





( 12 ) 



Thus the diffraction angle is approximately proportional to 
the input frequency. Figure 5 shows a classical optical ray 
trace diagram for the acoustooptical spectrum analyzer con- 
figuration. Figure 6 shows the corresponding optical plane 
schematic. In the output transform plane, position is 
approximately related to frequency by 



F = transform lens focal length. 

At this point, by considering the input signal as a 
composite of frequencies each of which will cause a deflected 
beam in accordance with the presented interaction theory, 
an approach to the construction of a model may be undertaken. 

Schulman [Ref. 18] and Vatz [Ref. 19] have proposed a 
set of optical block diagram symbols which show the process 
elements and indicate the mathematical function each element 
performs. Figure 7 shows the AOSA using these process 
elements . 

The laser light source S L , and beam expander lenses, T^ 
and T , are shown as providing an input to the Bragg cell 




(13) 



where 
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Figure 5. Classical Ray Trace Diagram for an Acoustooptical 
Spectrum Analyzer 
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Figure 6. Optical Plane Schematic for an Acoustooptical 
Spectrum Analyzer 
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Figure 7. Block Process Element Diagram for an Acoustooptical Spectrum Analyzer 
[after Refs. 18 and 19] 



with transverse light distribution a(x) of general form, an 
arbitrary profile modified by spatial noise, lens aberrations, 
and Bragg cell aperture. The input signal f(t) is converted 
to an acoustic wave in the E/0 block which represents the 
piezo-electric transducer. This acoustic traveling wave 
f(x-vt) propagates through the interaction region. 

In the Bragg cell interaction the light modulation can 
be considered a multiplication process in which the diffracted 
light power is proportional to the product of the light inten- 
sity and the acoustic wave power. 

i(x) = a(x) f(x-vt) (14) 

where 

x = v t (15) 

s 

The transform lens T produces the Fourier transform. 

I (s) = y {a (x) f(x-vt)} 

= A (s) * F ( s ) exp [ j 2 tt S (v t) ] (16) 

The transform is the convolution of the Fourier transform 
of the illumination function and the transform of the acoustic 
wave. The exponential factor is a Doppler frequency shift. 

As the detector is a square law device, the output from the 
AOSA is the square of the magnitude of the intensity function. 

| A ( s) * F ( s) | 2 (17) 

It has been shown that the output of the AOSA is the 
square of the magnitude of the intensity function. The in- 
tensity function may be considered as the convolution of the 
spatially transformed illumination function and the spatially 
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transformed acoustic wave. Hecht [Ref. 20] treats the illumi- 
nation function as a spatial weighting function composed of 
a number of factors including optical amplitude profile, 
aperture size, and acoustic attenuation. The AOSA frequency 
plane intensity distribution would thus be proportional to 
the squared amplitude of the Fourier transform of the weighting 
function with respect to the spatial frequency. 



x = finite period of signal analyzed 
S = fT = the normalized frequency 

Further refinement in the basic model is made by the 

consideration of the effects of transducer efficiency and 

diffraction efficiency. These factors modify the frequency 

components in the acoustic wave. 

C. AOSA MODEL CONSTRUCTION 
1. Irradiance Profiles 



the Bragg deflector is a complex decision. It would appear 
desirable to fill the entire deflector aperture with a uniform 
intensity light below the saturation level of the AOSA. The 
irradiance profile in this case would be that of a plane wave. 
A uniform intensity filling the entire aperture will 
maximize acoustooptic interaction and provides maximum 



F (s) 




(18) 



where u = - j] 



The choice of the optical beam profile incident on 
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incident light to be contributed to the deflected beam 
components. Abrupt windowing will truncate the beam profile 
and cause sidelobes to appear. These sidelobes, if they are 
of sufficient magnitude, may cause erroneous components to 
appear at the image plane or mask low intensity components. 

To reduce the sidelobe levels, a less discontinuous 
transition in the absorption-transmission regions is necessary. 
A suitable choice would be a Gaussian irradiance profile with 
a spherical phase wavefront. This irradiance profile common- 
ly occurs as a laser output and has well understood propaga- 
tion behavior. Smith includes a handy propagation program as 
part of his work [Ref. 3] . This program calculates Gaussian 
beam phase front radius and spot size for propagation through 
any number of thin lenses spaced at arbitrary distances. 

Since the Gaussian incident beam profile would not fill the 
entire aperture with uniform intensity and assuming continued 
operation below the AOSA saturation level, less light inten- 
sity would be available for diffraction into representative 
beams. This will result in the loss of sensitivity and the 
possibility of an inability to detect low amplitude frequency 
components in the presence of strong frequency components. 

Other beam irradiance profiles are possible and 
might be employed for their specific deflected beam charac- 
teristics. Figure 8 taken from Ref. 15 shows the deflected 
beam profiles for various incident profiles. Such curves 
are useful for comparison of expected sidelobe levels. 
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Field Amplitude 




Figure 8 . Bragg Deflection Beams for Various Incident 
Irradiance Profiles (after Ref. 15) . 
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Hecht's [Ref. 5] implementation of the Gaussian 



irradiance profile is of the form: 



wi(t f r,T) 



exp 





( 19 ) 



where 



T = D/2W q 

D = effective aperture 

W = Gaussian waist size 
o 

A significant design question in the utilization of Gaussian 

irradiance profiles is the choice of the ratio of the aper- 
-2 

ture to the e intensity. This is discussed also by Smith 

[Ref. 3]. Beam spot size is defined as the perpendicular 

distance from the beam centerline to the point where the 

_2 

irradiance has dropped by a factor of e . To use the full 
information capacity of the Bragg cell, the spot size must 
be enlarged to illuminate the complete cell aperture. In 
this situation the aperture will severly truncate the Gaussian 
beam resulting in a profile similar to that of a sine 
function. Varying the ratio of the spot size to the aperture 
will alter the level of the sidelobes of the profile. Smith 
limited his spot size to avoid sidelobes and thus isolate 
their effect in his study of optical excision. 

Hecht expresses his spatial weighting functions in 
equivalent time variables. The linear relationship between 
the spatial and time variables and the analogous relationship 
between the spatial frequency and frequency provide great 
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flexibility in the construction of the computer program to 
model the AOSA. The Fourier and inverse Fourier transforms 
utilized in the AOSA program make no distinction as to the 
transforming variables. 

Only plane and Gaussian irradiance profiles were 
included in the AOSA program. It would be a simple matter 
to incorporate other profiles such as the Lorentzian or 
triangle beams. However, the frequency with which the more 
exotic profiles would be encountered does not merit the 
addition. Figure 9 shows the comparative image plane 
profiles for both irradiance profiles with a monotonic 70 MHz 
input signal. In Figure 9 the plane wave illumination chosen 
completely illuminated the aperture of the Bragg cell. The 
Gaussian illumination used an arbitrarily chosen 3.8 milli- 
meter spot size. Comparison shows a near 50% reduction in 
the level of the Gaussian sidelobes to the level of the plane 
wave sidelobes but with a tradeoff of a decrease in the peak 
intensity response. The units for the Image Plane Intensity 
axis in this figure and all succeeding AOSA program figures 
are arbitrary. 

Figure 10 shows the image plane intensity for an 
incident Gaussian beam with a spot size of 1.945 millimeters. 
This is the "optimum" spot size according to Smith's 
criteria for sidelobe suppression. The sidelobes are barely 
visible but again a tradeoff results in an additional decrease 
in peak intensity. Use of a Guassian beam with this spot 
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Figure 10. Image Plane Intensity Profile (Magnitude) 

Produced by a Monotonic 70 MHz Signal with an 
Incident Gaussian Beam and a Simple Transducer 
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size narrows the portion of the Bragg cell aperture receiving 
full illumination. The portion receiving the full illumina- 
tion is termed the effective aperture. Its effect will be 
discussed in a later section. Figure 11 shows for comparison 
the image plane intensity profile of a fully illuminating 
plane wave. 

2 . Aperture Size 

The universal feature of realistic weighting functions 
is that of a finite aperture corresponding to a finite period 
of signal t represented by the unity amplitude window function 
[Ref. 5] . 




where 

T = D/V s 

Early approaches to acoustooptical spectrum analysis modeling 
considered the input signal as a continuous wave. The effec- 
tive aperture of the selected Bragg deflector truncated the 
incident light, such that only the segment of acoustic wave 
within the optical aperture region was intercepted. That 
segment of acoustic wave corresponds to the finite period of 
signal being analyzed. 

Improvements in Bragg deflector manufacturing tech- 
niques have yielded high quality, low acoustic velocity 
deflector crystals. These crystals are of increasingly larger 
dimensions. The Rayleigh limit resolution, N, is defined as 
the total deflection angle divided by the angular spot size. 
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Figure 11. Image Plane Intensity Profile (Magnitude) Produced 
by a Monotonic 70 MHz Signal with an Incident Plane 
Wave and a Simple Transducer 
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The relation simplifies to 

N = ti x if (21) 

where 

ti = acoustic transit time across the optical aperture 

Af = deflector bandwidth 

For the AOD-70 deflector (assuming that full aperture illumi- 
nation is provided) 

ti = 2.56 x D (microsec) (22) 

where D is the optical aperture in centimeters. 

Using the AOD-70 parameters of a 40 MHz bandwidth and a maxi- 
mum t x of 10 microsec based on the full available aperture 
of 3.9 cm, calculation yields a Rayleigh resolution capacity 
of 400 spots. To fill the full aperture anamorphic optics 
such as cylindrical lenses or telescopes must be used. Using 
standard procedure for the detector array utilization each 
spot is focused to fall across two adjacent photosites. The 
Rayleigh resolution capacity of 400 spots provides for a 
maximum frequency resolution of .1 MHz for the full 40 MHz 
bandwidth. 

Using a Gaussian beam with spot size of 1.945 milli- 
meters to fully suppress the sidelobes in the image plane 
intensity profile results in an additional penalty in a 
decrease in the Rayleigh resolution capacity of the "time- 
bandwidth product , " as given by Equation 21 . The interaction 
length is no longer the Bragg cell aperture but rather the 
effective aperture, in this case estimated as the region 
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twice the spot size symmetrical about the cell spatial midpoint. 
This results in a smaller optical aperture, a smaller acoustic 
transit time, and a resultant decrease in the Rayleigh resolu- 
tion. Thus, while a choice of spot size of this magnitude 
will suppress sidelobes, it requires a possibly unacceptable 
tradeoff in frequency resolution. 

The larger apertures bring an additional consideration 
to the AOSA model construction. A typical application input 
signal pulse might be of one microsecond in duration while 
the total acoustic transit time for the AOD-70 is ten micro- 
seconds. A one microsecond pulse input to the deflector will 
persist in terms of the frequency spectra it contains and 
therefore, in the deflected beam components incident on the 
image plane for the full ten microsecond travel time. This 
effect of "time spreading" of the r-f signal is significant 
for its enhancement of the detection probability for short 
duration signals. In the application of AOSA techniques to 
Electronic Support Measures (ESM) receivers, this is of prime 
import. A 0.1 microsecond r-f signal, which is currently 
beyond the capability of many traditional technology base 
receivers, would be stretched to a 1.0 microsecond pulse 
greatly improving its detection probability. It should be 
noted that while this "time spreading" increases detection 
probability, it reduces the capability for accurate direct 
pulse width measurement, also a highly desirable factor in 
ESM application. 
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i 



In summary the aperture size weighting function 



imposes three constraints on the AOSA model. The physical 
parameter of the optical aperture is modeled by a simple 
unity amplitude window function. The Rayleigh resolution 
and thus frequency resolution is simulated by the technique 
of zero padding while applying the Fourier transforms to 
achieve a digital record length which will yield the appro- 
priate frequency resoltuion. Zero padding is a digital 
implementation technique used in Fourier analysis whereby 
the desired sampled signal data is entered followed by a 
data string of zeroes. This increases the record length to 
provide a desired resolution in the transformed domain with- 
out affecting the sampled data validity. The transformed 
domain resolution is simply the inverse of the record length. 
The finite signal length is modeled with a unity amplitude 
window function implicit in the input signal as it is digitally 
implemented. Since the period of signal analyzed may be 
shorter than the acoustic travel time of the optical aperture, 
it is desirable to analyze the deflection caused by inter- 
action of the pulse acoustic segment at various points as it 
traverses the optical aperture. A shiftable digital array 
was utilized controlled by a time delay variable; zero delay 
corresponding to a symmetric straddling of the acoustic pulse 
segment about the spatial midpoint of the optical aperture. 

This time delay variable forms the basis for a dynamic AOSA 
program to be discussed in a later section. 
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3. Acoustic Attenuation 



As the acoustic wave travels the length of the Bragg 
deflector, it is subject to acoustic attenuation. The ampli- 
tude of the acoustic wave decreases, with a resultant decrease 
in the ability to interact with the incident light wave. 

Hecht [Ref. 5] treats the acoustic attenuation as a spatial 
weighting function of the form; 

w(t,a,r) = exp (-at) (23) 

where a is a frequency dependent loss factor 



with 



a - 



a o f 



m 



20 log 1() e 



(24) 



• U.JJ 

a in ^ TT m 

o ysec GH 111 
z 

Values for a Q have been tabulated for common acoustooptic 

modes. Acoustic attenuation typically follows power law 

frequency dependence with exponent, m, near 2. 

The AOD-70 Bragg deflector is manufactured with high 

quality flint glass as the interaction medium and Indium 

bonded Lithium Niobate for the piezo-electric transducers. 

The a Q values for this type of deflector range from .1 to .15 

2 

(dB/microsecond GHZ ) . A plot of the more stringent condi- 
tion is shown as Figure 12. Calculation of the acoustic 
attenuation for the AOD-70 device parameters in terms of 
a normalized transmission factor was carried out. The 
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figure 12. Attenuation Factor as a Function of Frequency for Lithium Niobate 
Transducer Deflector 



results are plotted as Figure 13. These results show a near 
unity transmission for the region of concern, from 50 to 
90 MHz. Since the transmission loss was negligible for this 
application, an acoustic attenuation factor of Hecht ' s form 
was not implemented in the AOSA programs. 



higher frequency ranges, acoustic attenuation will become a 
significant factor. Figure 13 shows increased loss as the 
frequency nears 1 GHZ. Further discussion of the effect of 
acoustic attenuation on weak truncations (Gaussian beam) and 
strong truncations (uniform optical beam) is given by Hecht 
[Ref. 21]. 

4 . Diffraction Efficiency 

For a simple transducer the Bragg angle cannot be 
exactly matched by the incident collimated light beam for 
all frequencies over the bandwidth of concern. As a result, 
for Bragg regime deflection the intensity of the diffracted 
beam is a function of the applied frequency. The ratio of 
the diffracted beam intensity to the intensity of the incident 
beam is defined as the diffraction efficiency. An analysis 
of Bragg regime deflection was conducted by Reagan in Ref. 2 
and is summarized here. 

Hecht [Ref. 5] expresses diffraction efficiency in 
a simplified form as: 



If the use of the AOSA programs is envisioned in 




( 25 ) 
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Figure 13. Normalized Transmission as a Function of Frequency for Lithium Niobate 
Transducer Deflector 



with 








P L 
a 

4 cos 2 8 

o 



(26) 



and 

L = nA 2 cos 0 /A^ (27) 

o o o o 

where: Ii = diffracted beam intensity 

I = undeflected beam intensity 
A q = free space laser wavelength 
n = refractive index 
p = elastooptic coefficient 
p = mass density 
A q = acoustic wavelength 
v = acoustic wave velocity 
P = acoustic power 

cl 

H = acoustic beam height 

0 Q = Bragg deflection angle 

L = interaction length 

L q = characteristic interaction length 

F = normalized center frequency = 1 

= input frequency value normalized to the center 
frequency value f 

Equation 26 shows the quantity -n to be constant with 
respect to the input signal frequency. Thus the initial sin 
squared term of Equation 25 may be ignored when considering 
the bandshape of a normal acoustooptic beam deflector. A plot 
of Equation 25 for various L/L q ratios is shown as Figure 14. 
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Figure 14. Diffraction Efficiency Bandshape (after Ref. 2) 
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The 3dB bandwidth is given by Equation 28 and is also from 

Ref. 5 . 

Af = 1.8 f L /L (28) 

o o 

For application of a Bragg regime deflector as an 
acoustooptical spectrum analyzer, it is desirable to maintain 
a bandshape as near unity and as uniform as possible over 
the bandwidth of concern. The characteristic interaction 
length is shown by Equation 27 to be constant with respect to 
dimensions of the Bragg deflector. The trend for AOSA 
applied Bragg deflectors is that of increased size with an 
attendant increased interaction length L. This reflects 
both improvement in manufacturing techniques permitting the 
production of high quality large deflectors and the desire 
for increased resolution made possible with larger inter- 
action length. The increased interaction length, however, 
results in higher value L/L Q ratios with the observed 
effect of narrowed bandwidths. 

Attempting to increase bandwidth by shortening the 
interaction length is of limited benefit since the diffraction 
shifts from the Bragg regime to the Raman-Nath regime. 
Raman-Nath regime diffraction is unsuitable for this appli- 
cation as it supports many higher order diffraction modes 
resulting in interference created by additional beams. 
Increasing the refractive index or elastooptic coefficient 
of the interaction medium, decreasing the acoustic beam 
height, or decreasing the acoustic velocity or the mass 
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density of medium can improve the diffraction efficiency to 
a limited extent. An increase of acoustic beam power will 
increase diffraction efficiency. The increased power also 
tends to heat the interaction medium, creating temperature 
gradients which degrade resolution and information capacity. 
Increased difficulty in internal reflection damping is 
experienced, permitting unwanted signal persistence. 

The only effective manner in which to increase the 
bandwidth is to increase the cell's center frequency f Q . 

The maximum frequency which the input transducer can effec- 
tively couple to the medium is determined primarily by the 
bonding technique used to attach the piezo-electric strips 
to the interaction medium. Thus the ability to increase the 
center frequency is dependent on the input signal transducer 
Normally Bragg deflectors are limited to operation over one 
octave . 

In summary, for a simple transducer little can be 
done to improve diffraction efficiency or increase bandwidth 
Beam steering, however, can significantly improve both 
diffraction efficiency and bandwidth. Commonly employed 
beam steering techniques include first order planar and 
stepped arrays. Figure 15 shows a typical first order 
stepped array. Pinnow [Ref. 22] discussed the relative 
merits of simple, planar and stepped array transducers. 

The larger dimensioned Bragg cells deflectors are almost 
always of the stepped array type. The AOD-70 deflector is 
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Figure 15. First-order Stepped Array Transducer 
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a five element stepped array. Figure 16 is the manufacturer's 
supplied efficiency versus frequency response curve. 

The AOSA program has both a simple transducer and 
a stepped array approximation option. The simple transducer 
response is computed in accordance with isotropic theory 
based on the AOD-70 device full dimensions. This permits a 
comparison with the actual stepped array response. Figure 17 
shows the computed diffraction efficiency curves for various 
L/L q ratios. These were obtained by sweeping a single varia- 
ble frequency over the band of concern and extracting the 
magnitude value for the instantaneous frequency applied. In 
comparison with Figure 14, which contains the bandshapes for 
the diffraction efficiency theory, a rippling is noticeable 
at the lower frequencies. This rippling is an undesirable 
artifact of the digital Fourier transform application. The 
ripple amplitude is comparatively small and falls outside of 
region of concern; hence it was not corrected. 

The stepped array bandshape that was included in the 
model is an approximation based on the near uniform response 
that was experimentally measured for the AOD-70 device. A 
slight roll-off in response occurs at the endpoints of the 
band, but for a first order approximation a uniform window 
function utilizing the mean efficiency was employed. 

To compare the response of the models of the two 
transducers, an input signal consisting of an in-phase 
frequency comb with 5 MHz spacing was applied. Figure 18 shows 
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Figure 18. Image Plane Intensity Profile (Magnitude) Produced 
by a 5 MHz Incremental Comb Signal with an Incident 
Plane Wave and a Simple Transducer 
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the simple transducer response. Figure 19 shows the stepped 
array response. Comparison shows severe attenuation of the 
frequency components beyond 80 MHz and a lesser degree of 
attenuation from 50 to 60 MHz for the simple transducer. 

On the other hand the stepped array transducer shows a very 
uniform response but with a tradeoff in decreased overall 
intensity magnitude. 

5 . Transducer Efficiency 

The measure of the transducer's ability to convert 
the r-f input signal components to representative acoustic 
waves is termed the transducer efficiency. The transducer 
efficiency is dependent upon the transducer material, the 
bonding technique used to fasten the transducer to the inter- 
action medium, and the transducer impedance matching. A funda- 
mental tradeoff relationship exists between the flatness of 
the frequency response and the amount of power delivered by 
the transducer within the desired band. The power levels 
used to drive the transducer are relatively low, thus we can 
afford a power conversion inefficiency and seek as uniform 
a response as possible. All transducers operate at less than 
unity efficiency; however, for purposes of the model, it was 
assumed that the desired uniform response was of paramount 
importance. An ideal case would provide a transducer efficiency 
response with little power attenuation and a uniform amplitude 
for all frequencies within the range of the system application. 
Current technology has achieved near uniform response bandshapes 
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Figure 19. Image Plane Intensity Profile (Magnitude) Produced 
by a 5 MHz Incremental Comb Signal with an Incident 
Plane Wave and a Stepped Array Transducer 
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over bandwidths of 70 MHz, as is present in the AOD-70 cell. 
Chang [Ref. 23] discusses the state of current transducer 
technology . 

Based on the above the AOSA model was implemented 
with a unity amplitude window for the transducer efficiency 
function. The program structure permits easy input of non- 
unity bandshapes should this be desired. Extreme broadband 
applications of acoustooptical spectrum analysis will result 
in nonuniform frequency responses requiring such modifications 
before using the AOSA program. 

D. AOSA MODEL UTILIZATION AND EXTENSIONS 

The AOSA program was designed to model direct detection 
acoustooptical spectrum analyzer systems in general. In its 
current configuration it utilizes the component parameters 
of the NPS laboratory system with the addition of the AOD-70 
Bragg cell. The coding for the Fortran program was selected 
to maintain as much generality as possible, so as to permit 
application for analysis under other choices of components. 

An attempt was made to minimize the number of preliminary 
calculations a prospective user would have to perform prior 
to utilization of the program. For example, the program 
computes the Fourier transform lens focal distance given the 
deflector physical parameters, the laser characteristics, 
and the desired detector array spot size. A list of the 
required input physical parameters is given in Table II. 



; 
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TABLE II 



AOSA Modeling Parameters 



DCELL 

LW 

VS 

REFIND 

LAMBDA 

BW 

SPOT 

PULSE 

TDELAY 



effective optical aperture of Bragg cell 
Bragg cell interaction length 
acoustic velocity of interaction medium 
refractive index of interaction medium 
laser wavelength 

beam width of collimated Gaussian beam (3 dB) at 

the Bragg cell aperture 

desired detector array spot size 

r-f signal pulse duration 

time delay for examining signal pulse at varied 
points during the acoustic transit period 
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The remaining program variables are primarily concerned with 
control of the program including output formats and the var- 
ious subroutines utilized within the program. 

The main exception to the above is the created variable 
TDELAY , for time delay, which permits the user to examine 
the representative acoustic pulse at frozen intervals of time 
during its transit of the optical interaction aperture. This 
suggested the possibility of adding another dimension to the 
AOSA program, to create in effect a dynamic AOSA model which 
would examine the response of the acoustooptic spectrum 
analyzer periodically by looping through the basic static 
algorithm with decreasing time delays. The results are 
presented in a three-dimensional isometric or perspective 
projection. This permits the user to simulate the response 
of the image plane intensity as the acoustic pulse transits 
segments of the optical interaction region rather than at a 
single fixed point of time. To simplify the graphic output 
the present configuration of the dynamic AOSA program (called 
TAOSA) is set to examine the first fraction of a microsecond 
as the acoustic pulse enters the optical region. Computer 
simulations were conducted over the entire transit period 
and to verify the symmetry of response as the acoustic pulse 
exited the optical region. 

Both the static and dynamic programs provide the user 
with options to simulate various system configurations. 

Types of transducers, incident beam profiles, and output in 
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either relative magnitude or decibel form is available. The 
static program offers a choice of several levels of graphical 
output, from a coarse line plot to two Versatec offline 
plotting options. Currently both programs utilize a deter- 
ministic form for the input r-f signal. The programs could 
be easily modified to accept sampled r-f data. 

The following series of figures are typical of the outputs 
of the AOSA and TAOSA programs. They were selected to demon- 
strate the ability of the programs to model a variety of 
complexity in input signals, to demonstrate the quality of 
resolution achievable by an acoustooptical spectrum analyzer 
and for comparison with the known spectral content of the 
inputed signals. 

Figure 20 is the input r-f signal magnitude for a 70 MHz 
monotone. Figures 21 and 22 are outputs using the decibel 
option of the AOSA program. Figure 21 shows the relative 
sidelobe levels by interaction of the monotone with an inci- 
dent plane wave. For comparison Figure 22 shows the decrease 
in the relative sidelobe level produced by the monotone with 
an incident Gaussian beam of 1.045 millimeter spot size. 

Figures 23 and 24 show the dynamic response of the same 
monotone with the two wave types. The x axis on all TAOSA 
program graphics represents time as the acoustic pulse transits 
the Bragg cell aperture region. Time zero corresponds to the 
instant when the leading edge of the acoustic wavetrain 
crosses into the aperture region. Figure 23 shows the dynamic 
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Figure 20. Input RF Signal Magnitude for a 70 MHz Monotone 
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Figure 21. Image Plane Intensity Profile (Decibel) Produced 
by a Monotonic 70 MHz Signal with an Incident 
Gaussian Beam and a Stepped Array Transducer 
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Figure 22. Image Plane Intensity Profile (Decibel) Produced 
by a Monotonic 70 MHz Signal with an Incident 
Plane Wave and a Stepped Array Transducer 
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Figure 23. Image Plane Intensity Profile (Magnitude) Produced 
by a Monotonic 70 MHz Signal with an Incident Plane 
Wave and a Stepped Array Transducer (Dynamic Model) 
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response of the monotone with an incident plane wave. 

A much slower response is seen in Figure 24 showing the 
response produced by an incident Gaussian wave of 1.945 
millimeter spot size. The Gaussian response is shown to 
occur later in the cell aperture transit period and to last 
for a relatively short time when compared to the plane wave 
response. The symmetric roll-off is already visible at the 
rightmost edge of the figure. 

Figure 25 is the input r-f signal magnitude for a simple 
sinusoidally modulated FM signal with 3=8.0 and f = 5 MHz. 

3 is the modulation index and is defined as 

3 = (29) 

m 

Figure 26 shows the image plane intensity profile produced 
by the interaction of the signal shown in Figure 25 with an 
incident plane wave. Figure 27 shows the image plane inten- 
sity profile produced by the same signal interacting with an 
incident Gaussian beam of 1.945 millimeter spot size. 

Figures 28 and 29 are the analagous dynamic program outputs. 
Again the Gaussian beam interaction shows a slower and 
smaller response. The advantage of sidelobe suppression, 
however, is apparent in the extremely clear presentation of 
the FM intermodulation products. 

Figure 30 shows the image plane intensity profile produced 
by a linear chirp with an incident plane wave. Figure 31 
shows the profile produced by the interaction with a Gaussian 
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Figure 24. Image Plane Intensity Profile (Magnitude) Produced 
by a Monotonic 70 MHz Signal with an Incident 
Gaussian Beam and a Stepped Array Transducer 
(Dynamic Model) 
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Input RF Magnitude for an FM Modulation Signal 

with 3 = 8.0, f = 5 MHz 
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Figure 26. Image Plane Intensity Profile (Magnitude) Produced 
by an FM Modulation Signal, 6=8.0, f = 5 MHz, 
with an Incident Plane Wave and a Stepped Array 
Transducer 
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Figure 27. Image Plane Intensity Profile (Magnitude) Produced 
by an FM Modulation Signal, 3=8.0, f = 5 MHz, 
with an Incident Gaussian Beam and a Stepped Array 
Transducer 
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X axis - Time (ysec) 

Y axis - Frequency (MHz) 
Z axis - Magnitude 

(arbitrary) 



Figure 28. Image Plane Intensity Profile (Magnitude) Produced 
by an FM Modulation Signal, B = 8.0, f = 5 MHz, 
with an Incident Plane Wave and a Stepped Array 
Transducer (Dynamic Model) 
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Figure 29. Image Plane Intensity Profile (Magnitude) Produced 
by an FM Modulation Signal, 6= 8.0, f =5 MHz, 
with an Incident Gaussian Beam and a Stepped Array 
Transducer (Dynamic Model) 
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Image Plane Intensity Profile (Magnitude) Produced 
by a Linear Chirp with Af =30 MHz, f =50 MHz, 
t = 2 usee with an Incident Plane Wave°and a 
Stepped Array Transducer 
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Image Plane Intensity Profile (Magnitude) Produced 
by a Linear Chirp with Af =80 MHz , f =50 MHz, 
t = 2 ysec with an Incident Gaussian B§am and a 
Stepped Array Transducer 
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80 MHz, 



beam. For both the chirp parameters are Af = 
f = 50 MHz, and t = 2 microseconds. Both show high and 
low frequency ripples. The high frequency ripple is an arti- 
fact of the digital implementation of the linear chirp. The 
lower frequency ripple present is the effect of the funda- 
mental and sidelobe interaction. The Gaussian interaction 
shows the effects of its inherent slower response and degraded 
frequency resolution in the attenuation of magnitude over the 
lower halfband frequency components. 

Figure 32 shows the dynamic response of the same linear 
chirp and an incident plane wave. The severity of the inter- 
action of the fundamental and sidelobe components, as addi- 
tional frequency components enter the Bragg cell, manifests 
itself in a dynamic transient response. Only after a con- 
siderable interval does the response settle into a steady 
state displaying the expected chirp spectral profile. For 
comparison Figure 33 shows the same linear chirp interacting 
with an incident Gaussian beam of 1.945 millimeter spot size. 
In this case the sidelobe suppression is a distinct advantage 
producing a smooth steady state response almost from the 
onset of the interaction. It must be kept in mind that the 
interaction will continue to begin much later than that of 
the plane wave and thus the plane wave will achieve steady 
state earlier than the Gaussian. 

The TAOSA program also has the decibel output option. 
Figure 34 shows the relative decibel image plane intensity 
of a monotone of 70 MHz interacting with a plane wave. 
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Figure 32. Image Plane Intensity Profile (Magnitude) Produced 
by a Linear Chirp with Af = 80 MHz, f = 50 MHz, 
t = 2 ysec, with an Incident Plane Wave and a 
Stepped Array Transducer (Dynamic Model) 
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Figure 33. Image Plane Intensity Profile (Magnitude) Produced 
by a Linear Chirp with Af = 80 MHz, f =50 MHz, 
t = 2 ysec with an Incident Gaussian B§am and a 
Stepped Array Transducer (Dynamic Model) 
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Figure 34. Image Plane Intensity Profile (Decibel) Produced 
by a Mono tonic 70 MHz Signal with an Incident 
Plane Wave and a Stepped Array Transducer (Dynamic 
Model) 
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III. SUMMARY AND CONCLUSIONS 



A. SUMMARY 

This thesis has presented an analysis of Bragg diffrac- 
tion theory and direct detection acoustooptical spectrum 
analysis modeling techniques. Based on the study of the 
various modeling theories a suitable approach for a first 
order simulation was selected and used to develop a Fortran 
program modeling the direct detection acoustooptical spectrum 
analysis process. Extensions of the basic theory have been 
added to include the factors of diffraction efficiency, 
transducer efficiency, irradiance profiles of the incident 
laser illumination, the aperture size of the Bragg cell, 
and acoustic attenuation of the generated acoustic wave. 

A second program dynamically models the AOSA process incor- 
porating a time variable to simulate the response of the AOSA 
as the acoustic wave transits the interaction region. The 
programs were used to model the NPS laboratory AOSA using 
the parameters from that system. 

B. AREAS FOR FURTHER STUDY 

Further effort related to this thesis is required in the 
following areas. 

1 . Nonlinearities 

The theoretical basis for the AOSA model program does 
not take into account the nonlinear effects of crossmodulation 
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(the compression of one signal due to another) or inter- 
modulation (the interaction of two or more signals) . While 
the restriction of the operation of the AOSA to the Bragg 
regime eliminates much of the effects of these nonlineari- 
ties, the overall effect on the replication efficiency is 
still at issue. Further analysis is needed to determine 
if these factors seriously degrade the model's accuracy. 
Hecht discusses this issue in Ref. 21. 

2 . Time Varying Signals 

The ability of the AOSA to respond to rapid time 
varying signals (which will depend on the speed of response 
of the detector array and the Bragg cell) needs to be 
examined. Currently the TAOSA program provides for a fixed 
delay for the simulation of overall system response time. 

3 . System Components 

The physical response and performance of the AOSA 
systems components warrant further study. Currently only 
a gross simulation of the system response time can be made 
by insertion of a fixed time delay. Closer examination of 
the speed of the detector array with particular emphasis on 
the uniformity of photosite response is needed. Integrated 
optical devices for AOSA and AOSE applications will use 
semi-conductor lasers. These lasers will need to be 
examined for the system effects produced by fluctuations in 
intensity and variations in coherence. The current HeNe 
laser used in the NPS system should be examined for its 
intensity and coherence fluctuations and resultant effects. 
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C. CONCLUSIONS 



Comparison of the AOSA model image plane intensity 
profiles with known AOSA system output spectral profiles 
demonstrates the AOSA program as a satisfactory first order 
simulation. The acoustooptical spectral analysis process 
yields high quality replication of the spectral content on 
input signals. Future coupling of the AOSA program with 
optical excision routines previously developed may provide 
additional insight into the application of acoustooptical 
spectral analysis and excision for use in wideband, high 
density signal environments. 
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APPENDIX A 



ACOUSTOOPTIC SPECTRUM ANALYZER (AOSA) PROGRAM 

This Fortran program calculates the relative Bragg 
deflected image plane intensity profiles for multifrequency 
input r-f signals to a modeled acoustooptic spectrum analyzer 
(AOSA) . The physical parameters used in the modeling of the 
AOSA are discussed within the program and are summarized in 
Table II of the thesis. 

The program provides several additional selectable 
variables : 

ITEST - provides for graphical and numerical or 
graphical only output 

IWAVE - provides for choice of plane or Gaussian 
incident light waves 

ITRANS - provides for choice of simple or first 

order stepped array transducer simulations 
IDB - provides for choice of graphical output in 
magnitude or relative dB formats 
I GRAF - provides for choice of line printer or 
Versatec graphical plotting 
IPLOT - provides for choice of PLOTG or DRAWP 
Versatec plotting subroutines 
Use of these variables is discussed within the program. 
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Additional alterable variables exist within the program 
related to the Fourier transform process. It is believed 
that the present configuration will suffice for most general 
applications . 

This program utilizes the PLOTG and DRAWP Versatec 
software subroutines, the line plot subroutine PLOTP, and 
FOURT, which calculates the forward and inverse Cooley-Tukey 
Fast Fourier Transform of multidimensional complex data. 
Prospective utilizers at the Naval Postgraduate School should 
refer to the W. R. Church Computer Center for detailed 
information. 
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ACOUSTO-OPTICAL SPECTRUM ANALYZER MODEL 
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APPENDIX B 



DYNAMIC ACOUSTOOPTIC SPECTRUM ANALYZER PROGRAM 

The dynamic acoustooptic spectrum analyzer program (TAOSA) 
is essentially a variation of the AOSA Fortran program. It 
computes the time delay to shift the acoustic pulse outside 
of the optical interaction aperture and then incrementally 
decreases this delay to simulate the action of the pulse as 
it transits the optical interaction region. The system 
response time may be simulated by adding a fixed time delay. 

The variables used in the program are discussed both 
within the program and in Appendix A. An additional variable, 
TINT, is added as an acceleration factor that varies the in- 
crement size based on the level of intensity of the incident 
laser beam. 

The results are presented as a three-dimensional isometric 
or perspective plot. Perspective users at the Naval Post- 
graduate School should refer to the W. R. Church Computer 
Center for detailed information on the three-dimensional 
plotting routine PLT3D1 in use. 
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AC GUSTO -OPTICAL SPECTRUM ANALYZER DYNAMIC MODEL 
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APPENDIX C 



ARBITRARY PROFILE SOLUTION DETAILS 
(after Ref. 15) 

The complex transmission coefficient T(0) is defined as 
T (0 ) = - a 1 ^ A ^ a ~ Det (cl) 

where 



Det = fa A - a B )fa A, 

y rrn * r 



2 aa 



' jfa A,, - a B.,J 
aa i ob 2 bo 



-a;a 2 (A , - B u )(k, - B, ) 

ab ah ba ba 



(C2 ) 



<*i " ITT" I' 31 * “ $ + - $> 2 + (Aftg^ 2 ] 1 ' 2 !. 

<*z * ,/3l ‘ “ ^ ~ $ )2 + {Mk ^.) 2 ] vz ), 

^- K h k S[ l+ 

x [exp (-ik^L) - expO'fcSjL)), 





X [exp(-iTtLL) - exp(i*iL>|, 



Au * fco\/Ticos0 » v / (2r/A) i €t - 3ij 

*?,-(* 3 < i -^-,) l/2 

*1.-(*fo-05> 1/2 



(C3) 

(C4) 



(C5 ) 

(C6) 

(C7) 

( C 8 ) 
(C9) 
(CIO) 
(Cll) 
(C12) 
(C13 ) 



M = index of Modulation of the permittivity of the Bragg 
cell medium 

£ ^relative permittivity of air ~1 
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e 2 = relative permittivity of the Bragg cell medium 



£ £ 

3 = 1 



9 = output angle 



k = 2tt/X 
o 

3 = (2tt/X) /z~ sin 9 (C14) 

o i 



8 = 3 - 2TT/d (C15 ) 

- \ o 
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